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The stability, electronic and optical properties of single-walled carbon nanotnbes (CNTs) and 
graphene noncovalently functionalized with free-base tetraphenylporphyrin (TPP) molecules is ad¬ 
dressed by density functional theory calculations, including corrections to dispersive interactions. 

We study the TPP physisorption on 42 CNT species, particularly those with chiral indices (n,m), 
where 5 < n < 12 and 0 < m < n. Our results show a quite strong tt-tt interaction between TPP 
and the CNT surface, with binding energies ranging from 1.1 to 1.8 eV, where higher energies can be 
associated with increasing CNT diameters. We also find that the TPP optical absorptions would not 
be affected by the CNT diameter or chirality. Results for the TPP physisorption on graphene show 
a remarkable stability with binding energy of 3.2 eV, inducing a small redshift on the 7r-stacked TPP 
absorption bands. The strong graphene-TPP interaction also induces a charge transfer from TPP 
to graphene, indicating a n-type doping mechanism without compromising the graphene structure. 


I. INTRODUCTION 

Noncovalent functionalization of single-walled carbon 
nanotubes (CNTs) with photoactive molecules is becom¬ 
ing a promising technique to explore functional materi¬ 
als for light-harvesting or optoelectronic applicationsii^— 
The coupling of the exceptional CNT transport proper¬ 
ties with the optical properties of functional dyes, like 
porphyrins with a strong absorption in the near infrared 
and visible, make of CNT/dye complexes good candi¬ 
dates for sensitive nanoscale devices with potential ap¬ 
plications in biomedical imaging^rJi and hybrid organic- 
inorganic photovoltaic devicesJ^ However, the strength 
of the TT-stacking interaction with respect to the CNT 
structural parameters, like the diameter and chirality, is 
still lacking. In addition, a better understanding of the 
CNT/dye optical properties considering changes in the 
dye structure due to the 7r-conjugation is needed. 

Recently, the kinetics and thermodynamics properties 
of non-covalently bound CNT/porphyrin oligomers have 
been investigated by UV/visible spectroscopy and fluo¬ 
rescence titrationJ^ It was reported that the affinity of 
the CNTs increases sharply with the porphyrin coverage, 
showing strongest binding energies for semiconducting 
CNTs, particularly those with chiral indices (7,5) and 
(8,6). In addition, others works have also investigated 
a supposed selective interaction of porphyrins towards 
semiconducting CNTs, suggesting that the semiconduct¬ 
ing and metallic CNTs would have significantly differ¬ 
ent surface properties These works have also sug¬ 
gested that noncovalent functionalization of a mixture 
of CNTs with porphyrins may be an effective method 
for the separation of semiconducting CNTs from metallic 
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CNTs. Previously, we reported physisorption properties 
of iron porphyrin on both metallic and semiconducting 
CNTs with similar diametersJ^ Our results showed a sur¬ 
prising strong TT-stacking interaction, but with negligible 
energy differences. Similar results were found for the ph¬ 
ysisorption strength of free-base and zinc porphyrins on 
the same semiconducting CNTfiid^ showing no clue on 
the suspected selectivity. 


Considering the increasing interest in the noncovalent 
complexation of CNT-based materials for optoelectronic 
applications, we carry out an in-depth investigation on 
the stability and optical response of free-base tetrapheny- 
porphyrin (TPP) molecules 7r-stacked on single-walled 
CNTs with different chiralities and diameters, and also 
on graphene. Our purpose is to give a theoretically in¬ 
sight into the binding strength and optical properties of 
these compounds as a function of the CNT structural 
parameters. Our results show that the diameter instead 
chirality would be the relevant CNT parameter to ex¬ 
plain the unusual strength of the TPP 7r-stacking inter¬ 
action. In addition, the TPP binding energy on metallic 
and semiconducting CNTs does not exhibit important 
differences that might suggest distinct interactions. Re¬ 
garding optical properties, we hnd that transition bands 
of the TT-stacked TPP remain at almost the same energy 
position than those found in the isolated TPP, being in¬ 
dependent of the CNT structural parameters. For the 
TPP physisorption on graphene, we find a very strong 
binding energy, of about 3.2 eV, providing an upper limit 
for the stability of CNT-TPP compounds. We also find a 
p-type doping on graphene induced by the TPP adsorp¬ 
tion. 
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II. THEORETICAL APPROACH 


Our density functional theory (DFT) calculations were 
performed using the SIESTA ab initio package;^ which 
employ norm-conserving pseudopotentials and localized 
atomic orbitals as basis set (double-^, singly polarized). 
The physisorption of TPP on both the CNT sidewall and 
graphene is assessed by van der Waals density functional 
(vdW-DF) as proposed by Dion et al^ This approach 
has been successfully applied to describe the dispersive 
interaction of aromatic molecules on the graphite sur¬ 
face, showing good agreement with available experimen¬ 
tal datai^ 

In recent experiments, the functionalization of CNTs 
with TPP molecules has been achieved by the micelle 
swelling methods in water suspension^i where the CNT- 
TPP compound stays in the micelle core. Because of the 
hydrophobic character of this core, it is not expected wa¬ 
ter molecules surrounding the compound, which can be 
considered close to a vacuum situation. In the present 
work we study the TPP physisorption 42 different CNT 
species, particularly those with chiral indices (n,m), with 
5 < n < 12 and 0 < m < n. The CNT-TPP 
compounds are studied within the supercell approach, 
with periodic boundary conditions along the nanotube 
axis {z direction). We use unit cells with a volume of 
(30 X 30 X £ao) , where oq is CNT lattice constant and i 
a factor to keep a minimum distance among TPP images 
along the CNTs, chosen to be of 10 A. This procedure 
results in supercell lengths (iao) between 22-57 A, which 
can contain up to 754 atoms. A grid cutoff of 100 Ry 
and the F point were used for the real-space and fc-space 
integration, respectively. For the optical calculations, we 
use a 1 X 1 X 31 fc-point mesh, for incident light polarized 
along the CNT axis. The accuracy of these parameters 
were tested considering larger grid cutoff (150 Ry) and 
fc-point mesh (1 x 1 x 10). Negligible variations in the 
total energies, optical spectra, and band structures were 
found, ensuring that our results are converged. For the 
TPP physisorption on graphene, we use a squared unit 
cell with 6x6 periodicity and the F point for the fc- 
space integration. Here, a vacuum region of 20 A among 
graphene an their images was considered. For the optical 
response we use a fc-point mesh of 31 x 31 x 1, for inci¬ 
dent light polarized parallel to the graphene plane (x and 
y directions). The TPP binding energy on the substrates 
is calculated by the energy difference between adsorbed 
and separated constituents, considering corrections due 
to the basis set superposition error. The complexes were 
fully relaxed by conjugate gradient minimization until 
the forces on the atoms were less than 0.05 eV /A. 

The optical response of the functionalized CNTs and 
graphene is obtained through first order time dependent 
perturbation theory, by calculating the imaginary part of 
the dielectric function (£ 2 ). £2 gives us a first approach 
for the optical absorption coefficient and it is calculated 



FIG. 1: Top and front views of a free-base tetrapheny- 
Iporphyrin (TPP) physisorbed on the (8,7) CNT in the equi¬ 
librium geometry. Blue and white bars represent nitrogens 
and hydrogens, respectively, whereas green and red are car¬ 
bons. 


according to the equation; 

ei{^') = A J dk^|e-(T,(k)|r|^„(k))|=5(Ec-E„-M. 

C,V 

( 1 ) 

Here A is a constant that depends on the cell sizes; d'c 
and are the occupied and empty Kohn-Sham orbitals, 
respectively. The delta function represents the conserva¬ 
tion of energy, which is described by a gaussian function 
with a smearing of 0.06 eV. 

It is important to note that DFT calculations fail in de¬ 
scribing electron-hole and electron-electron interactions, 
the so called many-body effects. These effects are respon¬ 
sible for the formation of excitons and the quasiparticle 
excitationj2^ Therefore, an accurate description of pho¬ 
toexcitations in CNTs and graphene needs theories be¬ 
yond DFT, like those based on the GW-Bethe-Salpeter 
equation (GW-BSE)i^ However, these calculations are 
not possible right now, considering the size of the sys¬ 
tems under study. Thus, £2 can give us a first approach 
for the optical response of 7r-stacked TPP on CNTs and 
graphene. The main difference of DFT calculations with 
respect to GW-BSE is a redshift in the absorption spec¬ 
trum. For instance, in graphene this redshift is found to 
be about 0.6 eV 
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III. RESULTS AND DISCUSSION 
A. Energetic and structural properties 

We first discuss different adsorption geometries for the 
TPP molecule on the CNTs. Figure [T] shows the equi¬ 
librium geometry of a TPP adsorbed on a (8,7) CNT. 
Here, the N-H bonds of TPP are oriented perpendicu¬ 
lar to the CNT axis, that we called position (i). Two 
other positions for TPP are likely to be found: with the 
N-H bonds parallel to the CNT axis (ii), and with the 
N-H bonds forming an angle of 45° with respect to the 
CNT axis (hi). Although position (i) is found to be the 
most stable, the other positions have total energies within 
0.3 eV, suggesting that all three are equally probables. In 
the equilibrium geometry, the TPP phenyl groups exhibit 
a small rotation, which come with a molecular twisting. 
To compare the stability of different CNT-TPP com¬ 
plexes, only position (i) will be considered throughout 
this work. The adsorption distance is measured between 
the CNT surface and the N atoms of TPP, as shown in 
Figure [H Table U lists TPP binding energies and bind¬ 
ing distances for all the complex under study. We also 
include the DFT bandgaps for the semiconducting com¬ 
plexes. Those with small bandgaps (~0.1 eV) are consis¬ 
tent with empirical model predictions of En as well as 
with experimental data<^ We observe that binding en¬ 
ergies and adsorption distances vary in ranges given by 
l.I < Eb < 1.8 eV and 3.0 < d < 3.3 A, respectively. We 
note that stronger TPP 7r-stacking interactions result in 
larger adsorption distance. This can be understood by 
look at the TPP phenyl groups, whose H atoms become 
closer to the CNT surface for large diameter CNTs, as 
can be seen in Figure[TJ To estimate an upper limit for Eb 
and d, we calculate the TPP physisorption on graphene, 
which can be considered as a CNT with an infinite diam¬ 
eter. We find Et, = 3.32 eV and d = 3.42 A. Although the 
binding distance is similar to those found on CNTs with 
D < 16 A, the binding energy is considerably larger, with 
almost twice the value found on the CNTs. This can be 
understood because of the larger graphene-TPP overlap 
area, which tend to maximize the 7r-stacking interaction. 

Figure [2] shows the CNT-TPP binding energy as a func¬ 
tion of the CNT diameter. We observe an almost linear 
increase in the TPP attachment strength with the CNT 
diameter, which in some way confirms our assumption 
that larger CNT-TPP overlap area tends to increase the 
binding energy of the molecule. However, this general 
tendency is not followed in some cases, as can be seen 
in Fig. [21 For instance, (i) the larger TPP binding en¬ 
ergy found for the (5,0) CNT (d? = 4 A) with respect 
to the (6,0) CNT (I) = 5 A), by about 0.2 eV, and (ii) 
the abrupt changes in the TPP binding energy profile go¬ 
ing from (6,0) to (6,1), (7,5) to (7,6), and (8,5) to (8,6). 
We attribute these results to the particular atomic ge¬ 
ometry of the CNT atoms just below the TPP molecule. 
We observe that the CNT-TPP binding distance tends 
to increase when the H atoms of the phenyl groups near- 


TABLE I: Adsorption distance (d), binding energy (Eb), and 
bandgap energy (Eg) of CNT-TPP complexes. D and 6c are 
the CNT diameter and chiral angle, respectively. Results for 
TPP on graphene (G-TPP) are also included. 


Complex 

d(A) 

^c(deg) 

EgieV) 

d(A) 

EbieV) 

(5,0)-TPP 

3.902 

0.0 

0.19 

3.000 

1.22 

(5,1)-TPP 

4.526 

8.9 

0.00 

3.052 

1.17 

(5,2)-TPP 

5.083 

16.1 

0.00 

3.106 

1.18 

(5,3)-TPP 

5.686 

21.8 

1.15 

3.124 

1.25 

(5,4)-TPP 

6.315 

26.3 

1.04 

3.142 

1.25 

(5,5)-TPP 

7.055 

30.0 

0.00 

3.064 

1.35 

(6,0)-TPP 

4.921 

0.0 

0.00 

3.117 

1.11 

(6,1)-TPP 

5.118 

7.6 

0.40 

3.084 

1.24 

(6,2)-TPP 

5.837 

13.9 

0.67 

3.116 

1.21 

(6,3)-TPP 

6.401 

19.1 

0.02 

3.180 

1.21 

(6,4)-TPP 

7.025 

23.4 

1.02 

3.079 

1.35 

(6,5)-TPP 

7.704 

27.0 

0.87 

3.161 

1.36 

(6,6)-TPP 

8.362 

30.0 

0.00 

3.114 

1.39 

(7,0)-TPP 

5.685 

0.0 

0.24 

3.091 

1.18 

(7,1)-TPP 

6.115 

6.6 

0.00 

3.056 

1.23 

(7,2)-TPP 

6.608 

12.2 

0.84 

3.134 

1.25 

(7,3)-TPP 

7.165 

17.0 

0.86 

3.183 

1.24 

(7,4)-TPP 

7.754 

21.1 

0.03 

3.166 

1.35 

(7,5)-TPP 

8.372 

24.5 

0.85 

3.138 

1.42 

(7,6)-TPP 

9.038 

27.5 

0.75 

3.142 

1.50 

(7,7)-TPP 

9.703 

30.0 

0.00 

3.182 

1.48 

(8,0)-TPP 

6.468 

0.0 

0.60 

3.159 

1.30 

(8,1)-TPP 

6.511 

5.8 

0.78 

3.106 

1.35 

(8,2)-TPP 

7.384 

10.9 

0.09 

3.143 

1.37 

(8,3)-TPP 

7.934 

15.3 

0.96 

3.211 

1.40 

(8,4)-TPP 

8.490 

19.1 

0.77 

3.163 

1.45 

(8,5)-TPP 

9.103 

22.4 

0.02 

3.192 

1.46 

(8,6)-TPP 

9.749 

25.3 

0.74 

3.298 

1.35 

(8,7)-TPP 

10.410 

27.8 

0.66 

3.197 

1.42 

(8,8)-TPP 

11.105 

30.0 

0.00 

3.110 

1.54 

(9,0)-TPP 

7.269 

0.0 

0.15 

3.213 

1.33 

(9,1)-TPP 

7.685 

5.2 

1.01 

3.145 

1.28 

(9,2)-TPP 

8.165 

9.8 

0.73 

3.248 

1.44 

(9,9)-TPP 

12.467 

30.0 

0.00 

3.205 

1.54 

(10,0)-TPP 

8.027 

0.0 

0.74 

3.179 

1.41 

(10,1)-TPP 

8.458 

4.7 

0.06 

3.100 

1.44 

(10,2)-TPP 

8.938 

8.9 

0.85 

3.126 

1.44 

(10,10)-TPP 

13.929 

30.0 

0.00 

3.222 

1.59 

(12,0)-TPP 

9.638 

0.0 

0.10 

3.180 

1.46 

(12,1)-TPP 

10.042 

4.0 

0.77 

3.163 

1.50 

(12,2)-TPP 

10.495 

7.6 

0.61 

3.123 

1.49 

(12,12)-TPP 

16.634 

30.0 

0.00 

3.180 

1.82 

G-TPP 

CXD 

- 

0.00 

3.424 

3.32 


est neighbors to the CNT surface, lie exactly above nan¬ 
otube C-atoms. This effect is due to the steric repulsion 
between them, which also explain the broad variation on 
the binding distances between the CNT surface and the 
TPP N-atoms, as shown in Figure |3l 

Figure|4]shows the CNT-TPP binding energy as a func¬ 
tion of the CNT chiral angle. Apparently, higher chiral 
angles would induce stronger TPP attachment. How¬ 
ever, if we compare the TPP attachment on CNTs with 
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FIG. 2: Binding energy of a TPP molecule physisorbed on 
different CNTs as a function of the CNT diameter (D). 



d{A) 


FIG. 3: Binding energy of a TPP molecule physisorbed on dif¬ 
ferent CNTs as a function of the CNT-TPP binding distance 
id). 


similar diameters and binding distances with those plot¬ 
ted as a function of the chiral angle, we observe that 
they have also similar binding energy, suggesting that the 
CNT diameter instead chirality is the relevant parameter 
involved in the strong TPP binding energy. For instance, 
if we take the CNTs with chiral indices (5,4) and (7,2), 
they show almost the same TPP binding energy when 
plotted with respect to both the CNT diameter (Fig. [2]) 
and adsorption distance (Fig. [3]). However, they also 
show the same binding energy when plotted with respect 
to the chiral angle (Fig. |T|). The same behavior can be 
verified for the pairs of CNTs (5,5)-(6,4) and (6,5)-(7,4). 



(deg) 


FIG. 4: Binding energy of a TPP molecule physisorbed on 
different CNTs as a function of the CNT chiral angle (Oc). 


These results suggest that the strength of the TPP in¬ 
teraction on the CNTs would not be related to the CNT 
surface geometry or chirality, but to the CNT diameter. 

The adsorption energy of benzene and naphthalene on 
CNTs (7-8 A diameter), using electron-correlated Moller- 
Plasset (MP2) method reports binding energies for ben¬ 
zene and naphthalene of 0.18 and 0.48 eV, respectively!^ 
This results suggest an overestimation of the TPP ad¬ 
sorption energy on CNTs. However, temperature pro¬ 
grammed desorption experiments performed in ultrahigh 
vacuum have measured the binding energy of polyaro¬ 
matic hydrocarbons (PAHs) adsorbed on graphite, show¬ 
ing quite strong desorption energiesIn that work, ben¬ 
zene, naphthalene, coronene, and ovalene adsorbed on 
graphite show desorption energies of 0.50, 0.85, 1.4, and 
2.1 eV, respectively. For the case of benzene and naph¬ 
thalene, the experimental values almost double the MP2- 
theory values. Even considering the CNT curvatures, the 
TT-stacking interaction appears to be much more larger 
than the usually expected. In addition, a recent study of 
PAHs adsorption on graphene shows that the vdW-DF 
calculations exhibit excellent agreement with the avail¬ 
able experimental data, while empirical and semiempir- 
ical calculations only preserve the correct trends In 
summary, the work concludes that the adsorption inter¬ 
action of TT-conjugate systems on graphene are a complex 
combination of dispersive and electrostatic interactions, 
which would explains to some extent our results for the 
large adsorption energies of TPP on both graphene and 
CNTs. 
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FIG. 5: Imaginary part of the dielectric function for semi¬ 
conducting CNTs with a vr-stacked TPP, for incident light 
polarized parallel to the tube axis. Red and black spectra 
indicate the optical response of CNT-TPP and the isolated 
CNT, respectively. The dashed lines indicate the position of 
Q and B bands of the isolated TPP. 


FIG. 6: Imaginary part of the dielectric function for semi¬ 
conducting CNTs with a vr-stacked TPP, for incident light 
polarized parallel to the tube axis. Red and black spectra 
indicate the optical response of CNT-TPP and the isolated 
CNT, respectively. The dashed lines indicate the position of 
Q and B bands of isolated TPP. 


B. Optical properties 

Noncovalent functionalization of a single-walled (6,5) 
CNT with TPP molecules have been reported in aqueous 
solution by mean of the micellar suspension methodi^ 
The optical absorption of (6,5)-TPP complexes shows 
resonances at 2.357 and 2.831 eV, which are associated 
to the Q and B (Soret) TPP bands, respectively;^ Our 
results for the optical response of the (6,5)-TPP complex 
for a polarization parallel to the tube axis, show clear 
transitions at 1.60 and 2.04 eV, which we associate to 
the TPP Q and B bands, respectively. The difference in 
energy between these peaks (0.44 eV) is almost the same 
than those observed in the experiments (0.474 eV)4 How¬ 
ever, the theoretical peaks show a redshift of about 0.8 eV 
with respect to the experimental ones, while the shape 
of the absorption spectra would be preserved. This shift 
is associated to many-body effects previously discussed. 
Therefore, through the £2 calculation it is possible to in¬ 
vestigate, qualitatively, changes in the optical properties 
of different CNT-TPP complexes in term of the CNT 
structural parameters, like diameter or chirality. 

Figure [S] and [S] show £2 as a function of the photon 
energy, for incident light polarized parallel to the tube 
axis, for the semiconducting CNT-TPP complexes. The 
red and black spectra indicate the optical response of 
CNT-TPP and the isolated CNT, respectively. The ver¬ 
tical dashed lines indicate the position of the TPP Q 
(1.60 eV) and B (2.04 eV) bands. The optical spectra 
clearly show the absorption of the 7r-stacked TPP at al¬ 
most the same position than the isolated TPP, indicat¬ 
ing that the TPP absorption properties would be pre¬ 
served in the complexes, being independent of the CNT 


diameter or chirality. However, we must note that the 
present results were obtained considering only one TPP 
molecule on the CNTs. Increasing TPP coverages may 
change this conclusion. The optical response of com¬ 
plexes considering metallic and semimetallic CNTs (not 
shown) exhibit the same properties. Another interesting 
feature is the resonance between the CNT E 22 transition 
with the B band for the CNTs: (5,4), (6,1), (7,3), and 
(10,0), and with the Q band for the CNTs: (6,1), (7,5), 
(7,6), and (7,2). This tuning between E 22 and the TPP 
bands would increase the emission of the corresponding 
complex. 

Figure [7] shows the electronic band structure and den¬ 
sity of states (DOS) of the graphene-TPP compound 
described in a squared reciprocal lattice. The vertical 
arrows show the optical transition between the frontier 
molecular orbitals of the 7r-stacked TPP, which corre¬ 
spond to the dispersionless subbands. Here, Qx {Qy) 
is the transition from HOMO to LUMO (LUMO-fl), 
while Bx (By) is the transitions from HOMO-1 to 
LUMO (LUMO-fl). These absorptions arise from tt — 
TT* transitions and can be characterized approximately 
by considering only the four molecular levels above 
mentioned In pristine graphene, the valence-band 

maximum and conduction-band minimum touch at one 
point in the k space, the Dirac point, which occurs at the 
Fermi energy. This indicates that graphene is a zero-gap 
semiconductor. In Fig. [3 we observe that the Fermi en¬ 
ergy (the dashed line) is shifted upward with respect to 
the Dirac point energy, which indicates a charge transfer 
from TPP to graphene, revealing a n-type doping mech¬ 
anism. Therefore, graphene would change its electronic 
characteristic induced by the strong TPP complexation. 
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As the TT-stacking molecule do not alter the graphene 
atomic structure, their transport properties would be 
preserved. Similar results have been reported for other 
aromatic molecules physisorbed on graphenei^ 



FIG. 7: Band structure and density of states of graphene with 
a TT-stacked TPP represented in a squared reciprocal lattice. 
The vertical arrows indicate the Q and B transitions between 
TT-stacked TPP frontier orbitals. The black DOS corresponds 
to pristine graphene. 

Figure[8]shows £2 of a TPP 7r-stacked on graphene, The 
G-TPP compound, considering an in-plain polarizations 
of incident light. For pristine graphene, we find a strong 
absorption in the infrared, below 1 eV, and a less intense 
absorption in the ultraviolet, at around 4 eV, which are 
attributed to transitions among tt and tt* bands. These 
results are in qualitative agreement with calculated ab¬ 
sorption spectra using the GW-BSE approach as well 
as with experiments;^!^ unless the redshift previously 
discussed. For functionalized graphene, we observe two 
intense absorptions at around 1.57 and 1.96 eV, corre¬ 
sponding to the Qxy and B^y bands of the 7r-stacked 
TPP. These peaks show a redshift with respect to those 
of the isolated TPP of 0.05 eV, which is associated to the 
distortion of the TPP phenyl groups due to the strong 
TT-stacking interaction. Additional peaks at around 2.6 
and 3.1 eV are also associated to the TPP molecule. 


IV. SUMMARY AND CONCLUSIONS 

We have theoretically investigated the stability and op¬ 
tical absorption properties of a free-base TPP molecule 
TT-stacked on single-walled GNTs, as well as on graphene. 
42 GNTs with chiral indices where 5 < n < 12 

and 0 < m < n, have been studied. In term of the di¬ 
ameter, they range within 3.9 < D < 16.6 A. The TPP 
physisorption is described by van der Waals density func¬ 
tional while the optical properties through the imaginary 



FIG. 8: Imaginary part of the dielectric function for graphene 
with a TT-stacked TPP, for in-plain polarization of the incident 
light. The dashed lines indicate the Q and B bands of the 
isolated TPP. 


part of the dielectric function in the linear optical re¬ 
sponse. 

We find a strong CNT-TPP Tr-stacking interaction, 
with binding energies and adsorption distance varying 
within 1.1 < Ef, < 1.8 eV and 3.0 < d < 3.3 A, re¬ 
spectively, where larger Ei, and d are obtained for in¬ 
creasing GNT diameters. Our results indicates that di¬ 
ameter instead chirality is the relevant parameter for 
the strong TPP physisorption, which would be origi¬ 
nated in the increasing overlap area between TPP and 
the GNTs. An upper limit for the GNT-TPP binding 
strength can be estimated through the TPP physisorp¬ 
tion on graphene, which can be viewed as an infinite- 
diameter GNT. We find very strong Tr-stacking interac¬ 
tions, with Eb = 3.32 eV and d = 3.42 A. However, this 
strong interaction is not enough to open an energy gap in 
the functionalized graphene, but it can induce a charge 
transfer from TPP, suggesting a n-type doping mecha¬ 
nism without compromising the graphene structure and 
transport properties. Concerning optical properties, we 
find that Q and B (Soret) TPP absorption bands remain 
at almost the same energy position after the physisorp¬ 
tion on the GNTs, being independent of the GNT struc¬ 
tural parameters. Whereas on graphene, the TPP bands 
show a redshift of about 0.05 eV, which are associated 
to the distortion of the phenyl groups due to the strong 
TT-stacking interaction. Finally, our results suggest that 
both compounds CNT-TPP and G-TPP exhibit remark¬ 
able stability while preserving the absorption properties 
of the chromophore, which may be of interest for light¬ 
harvesting and bio-labeling applications. 
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